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ABSTRACT 


A  simple  one -dimensional  theory  is  given  for  the  steady,  compressible, 
adiabatic  flow  of  a  perfect  gas  through  a  porous  plate.  The  Dupuit-Forchheimer 
relation,  valid  for  incompressible  flow,  is  replaced  by  an  isentropic  com¬ 
pression  when  the  gas  enters  the  plate  and  a  non-isentropic  sudden  enlarge¬ 
ment  process  when  it  exits.  A  generalized  form  of  Darcy's  equation  is 
used  that  is  applicable  to  adiabatic  flow.  It  retains  the  convective  term, 
which  is  necessary  if  the  flow  is  compressible.  An  important  consequence 
of  this  study  is  that  the  Mach  number  at  the  downstream  surface  may  be 
much  smaller  than  unity,  even  when  the  flow  through  the  plate  is  choked. 

As  the  pressure  ratio  across  the  plate  decreases,  the  flow  remains  choked, 
but  the  downstream  Mach  number  increases.  In  fact,  this  Mach  number  will 
be  greater  than  unity  for  a  sufficiently  small  pressure  ratio,  in  which  case 
the  downstream  flow  is  supersonic.  Thus,  a  wide  range  of  downstream 
Mach  numbers  trom  subsonic  to  supersonic  is  possible,  even  though  the 
flow  is  choked.  For  incompressible  flow,  the  volumetric  flow  rate  varies 
linearly  with  the  pressure  differential  across  the  plate.  The  equivalent 
compressible  relation  is  shown  to  consist  of  a  plot  of  upstream  Mach  number 
versus  the  pressure  ratio  across  the  plate  The  incompressible  result  can 
also  be  shown  on  this  plot;  it  differs  from  the  compressible  one,  except 
when  the  plate  is  thick. 
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NOMENCLATURE 


b, 

drag  coefficient; 

Bo* 

Darcy's  constant; 

d, 

average  pore  diameter; 

fi' 

functions  defined  by  Eqs.  (18)  and  (20); 

F(M2) 

fuactic  i  defined  by  Eq.  (14); 

G;{M2), 

Mach  number  functions; 

h, 

enthalpy; 

L, 

porous  plate  thickness; 

M, 

Mach  number; 

P» 

pressure; 

R, 

gas  constant; 

Rd, 

Reynolds  number: 

®, 

entropy; 

T, 

temperature; 

u, 

velocity; 

X, 

distance  in  the  direction  of  flow; 

X, 

drag; 

P. 

porosity; 

V* 

ratio  of  specific  heats; 

^1. 

viscosity; 

P. 

density. 
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INTRODUCTION 


For  the  vast  majority  of  problems  concerning  flow  o£  a  liquid  in  a 
porous  medium,  the  tiu-.d  is  considered  incompressible  (or  slightly  com¬ 
pressible)  and  the  convective  terms  are  neglected.  For  most  classes  of 
problems,  such  as  ^round  water  and  seepage  flow,  these  assumptions  are 
not  unreasonable,  and  have  been  used  in  investigations  of  this  type  [  i,  2] . 

When  the  fluid  is  a  compressible  gas  subject  to  large  pressure  gradients, 
such  assumptions  are  no  longer  valid,  and  a  different  methodology  must  be 
applied.  The  purpose  of  this-  paper  is  to  develop  a  simple  theory  for  .one¬ 
dimensional,  compressible  flow  of  a  gas  through  a  porous  plate. 

This  work  was  initially  motivated  by  recent  developments  in  high  speed 
reentry  technology.  In  an  attempt  to  simulate  ablation  off  the  surface  of  a 
reentry  vehicle,  Hartunian  and  Spencer  (  3]  introduced  the  concept  of  massive 
blowing  in  wind  tunnel  experiments  using  models  of  a  porous  material  such 
as  sintered  stainless  steel.  A  substantial  pressure  difference,  usually  greate 
than  100  to  1,.  is  established  across  the  wall  of  the  model.  A  large  mass  flow 
rate  through  the  wall  ensues  and  appreciable  density  gradients  occur. 

Out  interest  here  is  the  interaction  of  the  injected  gas  with  the  porous 
material,  not  the  external  flow  field  of  the  injected  and  free-stream  gases. 

In  order  to  concentrate  more  fully  on  the  physics  of  this  interaction,  the 
following  simplifying  assumptions  are  introduced: 

1.  Steady,  one -dimensional  flow  is  assumed. 

2.  The  physical  model  of  a  gas  flowing  in  a  straight,  frictionless 
duct  of  constant  cross  section  is  adopted.  Situated  in  the  duct 
is  a  porous  plate  of  uniform  thickness  (see  Fig.  1). 
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3,  A  perfect  g£S  is  assumed  with  a  constant  value  for  the  ratio  of 
specific  heats. 

4.  Cot.tinuum  mechanics  govern  thv,  'low. 

We  take  as  a  mathematical  model  for  the  viscous  flow  inside  the  plate 
the  formulation  applicable  to  flow  in  a  straight  duct  with  friction.  This 
description  is  intended  to  represent  average  flow  conditions  at  any  cross 
section  of  the  plate  and  does  not  imply  the  assumption  of  a  capillary  model 
in  which  the  pores  are  straight  ducts  transverse  to  the  plate.  The  theory 
associated  with  this  approach  is  well  established  and  widely  known.  (See, 
for  example,  Shapiro  [4],  whose  work  is  heavily  relied  upon  here. )  We 
follow  [4]  and  use  the  Mach  number  as  the  independent  variable,  which  yields 
a  formulation  that  is  physically  elegant  and  mathematically  simple.  An 
important  result  of  this  approach  is  that  momentum  considerations  can  be 
dealt  v/ith  separately,  as  has  been  done  in  the  subsequent  analysis. 

The  distinguishing  feature  of  this  analysis  is  the  generalization  of  the 
Dupuit-Forchheimer  relation  [  1,  2],  which  says  that  the  volumetric  flow  rate 
is  constant  across  the  surface  of  a  porous  medium.  The  fluid  thus  adjusts  to 
the  smaller  flow  area  available  to  it  in  the  medium.  In  incompressible  flow 
theory,  this  relation  is  used  regardless  of  whether  the  fluid  is  entering  or 
leaving  the  medium.  By  contrast,  instead  of  the  Dupuit-Forchheimer  relation, 
we  use  the  isentropic  assumption  when  the  fluid  is  entering  the  plate.  On 
exiting,  the  flow  is  assumed  to  adjust  to  the  larger  available  area  by  means 
of  an  irreversible  process. 


The  most  important  consequence  of  this  analysis  is  that  although  the 
flow  may  be  choked,  the  Mach  number  at  the  downstream  surface  is  generally 
not  unity.  (By  choking  we  mean  that  the  mass  flow  rate  through  the  porcus 
plate  is  a  maximum. )  Specifically,  as  the  pressure  ratio  across  the  plate 
decreases,  the  flow  becomes  choked,  but  the  downstream  Mach  number  may 
still  be  much  smaller  than  unity.  As  the  pressure  ratio  further  decreases, 
the  flow  remains  choked,  while  the  downstream  Mach  number  increases.  In 
fact,  this  Mach  number  will  become  greater  than  unity  for  sufficiently  small 
pressure  ratios.  In  other  words,  although  the  flow  through  the  plate  is  choked, 
a  wide  range  of  downstream  Mach  numbers  ranging  from  subsonic  to  super¬ 
sonic  is  nevertheless  possible.  In  addition,  a  limiting  pressure  ratio  and 
Mach  number  exist  beyond  which  solutions  are  no  longer  possible  according 
to  this  theory. 

Results  of  this  work  can  be  given  as  a  plot  of  upstream  Mach  rumber 
versus  the  pressure  ratio  across  the  plate.  This  is  equivalent,  for  the  com¬ 
pressible  case,  to  the  normal  incompressible  relation  of  the  volumetric  flow 
rate  versus  the  pressure  differential  across  the  plate.  The  incompressible 
result  can  also  be  shown  on  our  plot;  it  differs  from  the  compressible  one, 
except  when  all  Mach  numberr  are  small  or  when  the  plate  is  thick  (see 
below). 

We  formulate  the  problem  in  the  following  section.  The  case  when 
choking  does  not  occur  is  analyzed  first.  Momentum  considerations  are 
dealt  with  in  a  separate  section  and  the  validity  of  Darcy's  equation  for 


-3- 


bo*h  isothermal  and  adiabatic  flow  through  the  plate  is  discussed.  The 
subsequent  section  treats  the  choked-flow  case,  while  the  final  section 
discusses  a  number  of  important  aspects  of  the  theory,  including  the  gener¬ 
alization  of  the  plot  of  flow  rate  versus  pressure  differential. 


FORMULATION 


Consider  one -dimensional  flow  of  a  perfect  gas  through  a  plate  of 
uniform  porosity  P  (see  Fig.  i).  The  cross-sectional  area  of  the  plate  is 
taken  to  be  unity,  consequently,  the  area  available  to  the  flow  within  the  plate 
is  p,  which  must  be  less  than  unity.  A  streamtube  of  fluid  thus  contracts  in 
area  as  it  enters  the  plate  and  expands  in  area  as  it  leaves.  These  changes 
are  assumed  to  occur  between  locations  1  and  2,  and  3  and  4,  respectively. 
Both  area  changes  are  assumed  to  occur  in  a  short  distance  in  comparison  to 
the  plate's  thickness. 

Upstream  of  the  plate,  the  flow  is  at  a  low  subsonic  Mach  number, 
i.  e.  ,  «  1,  given  by 


m2  _  11 1 
1  " 


RT1(p1tt1)i 

_ 

YPl 


(1) 


where  u  is  velocity,  \  is  the  (constant)  ratio  of  specific  heats,  p  is  pressure, 
p  is  density,  P  is  the  gas  constant,  and  T  is  temperature.  Because  «  i, 
the  thermodynamic  quantities  in  equation  (1)  differ  negligibly  from  their 
stagnation  values.  Similarly,  in  all  ensuing  formulas  we  use  without  further 
statement  or  justification  the  approximation 

i  +  (Y  -  l)Mj/2  s  1 

Between  locations  1  and  2  we  assume  the  area  change  is  accomplished 
by  means  of  an  iseotropic  process;  the  flow  is  thus  analogous  to  that  in  the 
convergent  part  of  a  nozzle.  Associated  with  t*  e  flow  is  a  small  loss  of 
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stagnation  pressure,  but  this  loss  is  minute  compared  with  that  between 
locations  2  and  3  and  therefore  is  neglected. 

Between  locations  2  and  3  we  assume  the  flow  is  adiabatic  but  non- 
ieentropic  due  to  viscous -energy  dissipation.  *  Flow  conditions  are  described 
by  a  Fanno  curve,  which  combines  the  adiabatic  assumption  with  continuity. 
Such  a  curve  [4]  is  independent  not  only  of  the  form  of  the  drag  term  in  the 
momentum  equation,  but  of  the  momentum  equation  itself.  This  is  an 
i  nportant  point  to  bear  in  mind,  since  much  of  the  analysis  deals  with 
properties  of  the  Fanno  curve,  and  is  thus  valid  for  any  momentum  equation. 

On©  of  the  fundamental  properties  of  a  Fanno  curve  is  that  at  the  "nose" 
(see  Fig.  2)  the  entropy  is  a  maximum  and  the  Mach  number  is  unity. 
Associated  with  this  is  the  phenomenon  of  choking.  Specifically,  when  <  1, 
as  i3  the  case  here,  we  have  the  condition  M^  -  i,  which  plays  an  important 
role  in  the  subsequent  analysis. 

Between  3  and  4  we  assume  the  area  change  is  accomplished  by  a  sudden 
enlargement,  which  is  a  non-isentropic  process.  When  M^  <  1,  this  is  a 
compression,  i.  e.  ,  an  increasing  pressure,  which  results  in  a  lower  Mach 
number  at  4  than  at  3.  When  M^  is  near  unity  and  (5  is  small,  an  entropy 

*W«  could  also  assume  the  flow  between  2  and  3  to  be  isothermal.  However, 
[4]  (Chapter  6)  shows  that  both  assumptions  yield  similar  results  qualitatively 
and  quantitatively.  For  experiments  of  long  duration  involving  flow  through 
a  thin  plate,  we  adopt  the  adiabatic  assumption  as  the  more  realistic  one. 


increase  occurs  which  cannot  be  neglected.  Since  this  case  is  important,  we 
use  the  more  exact  irreversible  process  throughout  the  analysis,  even  though 
an  isentropic  process  would  otherwise  bs  feasible  when  <  1. 

These  considerations  are  conveniently  illustrated  by  a  Mollier  diagram 
(Fig.  2).  The  ordinate  is  the  enthalpy  h,  where  the  stagnation  enthalpy, 
is  constant  from  1  to  4,  since  the  entire  process  is  adiabatic.  Locations  i, 

2,  3,  and  4  are  state  points  on  the  diagram  and  hereafter  are  so  designated. 
States  1  and  £  ire  connected  by  a  Fanno  curve,  since  these  states  have 
the  same  mass  flux  rate  and  h^.  This  curve  is  to  the  right  of  the  curve  pass¬ 
ing  through  2  and  3  because  of  the  smaller  mass  flux  per  unit  area.  In  the 
P  =  1  limit,  both  curves  coincide;  as  p  decreases  from  unity  they  diverge. 

In  Fig.  2,  flow  conditions  between  i  and  2  are  given  by  an  isentropic  subsonic 
expansion,  flow  conditions  between  2  and  3  are  given  by  the  Fanno  curve,  while 
the  irreversible  subsonic  compression  occurs  between  3  and  4. 

Figure  2  and  the  process  between  states  3  and  4  are  valid  when  <  1 
and  in  the  limit  —  1.  When  =  i,  the  flow  chokes  and  the  analysis  must 
be  altered  to  allow  for  other  solutions  in  addition  to  the  -*  1  limit.  We 
shall  first  consider  the  <  1  case,  along  with  the  limit  i. 


M3  <  i  AND  M3  —  1  ANALYSIS 


States  1  and  4  are  related  by  the  Fanno  curve  [4] 


where 

GQ(M2)  =M2  [i  +  (y  -  1)  M2/sj. 

States  1  and  2  are  isentropically  related  by  [4] 

M2  =  p2  G^M2),  (2b) 


where 


/»» 2*  _  M 

G  (M  )  =  — - — — - - — - - - —  • 

r  ?  T(rH)/(Y-i> 

[i  +  (Y  -  1)  M^/2j 

States  3  and  4  are  related  by  the  momentum  equation  for  a  sudden  enlarge¬ 
ment  from  an  area  |3  to  unit  area.  This  equation  ie  (see  M.  problem  5.  23) 

(1  -  P)  p3  +  (P)(p3  +  p3  u2)  =  { i )  (p4  +  p4  u2),  (3) 
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which  is  valid  only  when  <  1,  or  in  the  limit  -*  1;  Bymeansof  the 
isentropic  relation  for  p.,/pj  and  the  jPanno  curve  relation.for  pV/p-ji 
equation  (3)  can  be  written  as 


The  three  equations  (2)  relate  four  Mach  numbers,-  and  consequently  an 

'  '  i*\.\  > ,  / 

additional  relation  is  necessary.  This  is  provided  by  conservafcidhdf  momen¬ 
tum  between  states  2  and  3,  which  is  given  later. 

t 

A  more  convenient  form  for  equation  (2c),  obtained  with  the  aid  of 
(2a)  and  (2b),  is 


M*  +  (y  -  l)M*/2  pfpMj)  +  (y  -  i|(pM|)  /2  ' 
(l  +  y  M^) 


A.  2 

[i  +  y  (^3)] 


(4) 
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We  deduce  from  this  that'M^  £  From  equation  (3)  we  have 


P4  14  V(PM|) 
p3  i  +  yM^ 


(5) 


and  consequently  p^  <  p^. 


Aiiicpnd.boiMxd-on  M^,  moreuseful  than  the  one  above,  is  obtained 
from  eqpation  {4)  and  the  condition  Mg  <  i,  and  is 


g2(m\) 


<  ^r1  (rhp) 


where 


(6) 


,  2  GMc) 

G2(M^)  *  - 

(l  + 


2  ‘ 


2  2  ' 

^Figure  3  is  a  plot  of  the  G.,(M  )  versus  M  for  y  =  (7/5).  Thus,  M4  is 

bounded,  where  this  bound,  M4 ,  is  obtained  from  (6)  by  replacing  the 

-~2 

inequality  by  an  equality  sign.  Figure  4  shows  M_,  versus  p  for  y  =  (1/5). 

2 

From  equation  (2a)  and  inequality  (6),  we  see  that  (pjMj/p4)  also 

has  an  upper  bound,  (p^Mj/p^r,  which  is  shown  as  a  function  of  M^in  Fig.  5 

~2 

for  y  =  (7/5).  For  a  given  value  of  p  there  is  a  maximum  value  (Fig.  4), 
which  in  turn  leads  to  a  maximum  value  (pjMj/p^/  (Fig.  5).  In  the  limit 
i,  we  have  =  M^  and  (p^Mj/p4)  =  (p^M^/p4).  When  Mg  =  1, 
however,  other  solutions,  which  are  discussed  later,  are  also  possible. 

It  is  clear  from  Fig.  4  that  M4  will  be  small  when  p  is  small.  Suppose 
both  Mj  and  M^  are  small.  This  does  not  imply  that  the  flow  is  incompressible, 
since  Mg  may  still  be  near  unity.  Figure  4  further  shows  (when  Mg  <  i)  that 
large  values  of  M^  cannot  be  attained  for  many  porous  media,  since  p  is  usually 
usually  less  than  0.  6. 


In  the  limit  when  all  Mach  numbers  are  small,  we  obtain  from 


equations  (2)  the  simple  result 


which  is  used  later. 


MOMENTUM 


foregoing  analysis  is  completely  independent  of  the  momentum 
between  states  2  and  3.  For  compressible  flow,  this 
written  as  ,(#ee  [4],  Chapter  8) 


dp  +  dX  =  -pu  du 


(8) 


The  drag  term  is 

dX  =~bpu2  dx  ,  (9) 

where  x  is  distance  in  the  direction  of  flow,  and  b  is  a  drag  coefficient  with 

-1 

dimension®  (length)  ,  considered  constant  when  any  integration  with  respect 
to  x  is  performed.  Equation  (9)  is  identical  to  that  used  in  [4]  (Chapter  6). 
with  b  s  (4f/D)v  and  was  chosen,  in  part,  for  this  reason.  A  comparison  with 
Darcy’s  law  (as  given,  for  instance,  in  [2],  p.  2)  shows  that  b  =  (2ti/puBq), 
whero  £>,»j  is  Darcy*®  constant  and  q  is  the  viscosity  of  the  gas.  A  more 
detailed  discussion  of  the  physical  implications  of  the  drag  coefficient  is  given 
in  the  last  section. 

The  Isothermal  Case 

Darcy's  equation  is  universally  used  as  the  momentum  equation  for 
incompressible  How  in  a  porous  medium.  Although  we  are  primarily  interested 
in  adibatic  flow,  it  is  useful  to  consider  the  isothermal  case  in  order  to  show 
the  relationship  of  Darcy's  equation  to  the  p work.  We  begin  by  pre¬ 
senting  the  usual  derivation  of  Darcy's  equation  for  isothermal  flow,  neglecting 


the  convective  term  pu  du.  We  then  present  a  more  precise  derivation, 
which  retains  the  convective  term  and  shows  that  Darcy’s  equation  is  not 
uniformly  valid  in  the  incompressible  limit.  Nevertheless,  it  is  also  shown 
that  Darcy's  equation  usually  differs  negligibly  from  the  correct  momentum 
equation.  As  we  shall  see  later,  the  same  type  of  nonuniformity  appears  in 
the  adiabatic  case. 

Consider  the  isothermal  case  neglecting  pu  du  and  the  area  changes  at 
the  plate's  surfaces.  The  area  changes  are  included  when  the  adiabatic  case 
is  analyzed.  We  therefore  obtain 

12  1  2 
p  dp  =  -  bppu  dx  =  -  bRT^pjU.)  dx 


and  by  integrating  between  1  and  4,  we  obtain  Darcy's  equation 


bL  = 


Pi  ’  P4  1  "  (P^Pt^ 

=  _ 


RTjtp^r 


yMj 


(10a) 


where  L  is  the  thickness  of  the  plate,  and  where  equation  (1)  is  used  to  obtain 

the  second  equality.  The  only  alteration  that  the  neglected  area  changes  make 
2 

is  to  introduce  a  (3  factor  on  the  right  side  of  equation  (10a) 

If  we  now  retain  the  pu  du  term,  and  still  assume  isothermal  flow,  we 
obtain  (see  [4],  Chapter  6) 


bL  =- 


73 


yM 


Y 


+  2  lr 


M 


t  V 


By  means  of  an  equation  for  similar  to  (1)*  this  becomes 


bL  - 


*  -  {P4/Pt> 


YM 


1 (P4/Pj[) 


(10b) 


which  differs  from  (10a)  by  the  logarithm  term.  Because  of  this  term,  Darcy's 

equation  is  not  a  uniformly  valid  approximation  in  the  incompressible  limit. 

(In  this  limit,  we  replace  by  equation  (1). )  Figure  6  shows,  however, 

that  the  absolute  value  of  the  logarithm  term  is  small  compared  to  the  value 

2 

of  the  right  side  of  (10a),  except  for  exceedingly  small  values  of  (P4/Pj)  , 
when  compressibility  should  not  be  neglected  in  any  case. 

The  Adiabatic  Case 

Next,  consider  adiabatic  flow  from  states  2  to  3,  and  in  parallel  to  the 
preceding  case,  examine  first  equation  (8)  without  the  pu  du  term.  This 
equation  can  be  written  as 


b  dx 


•  iljJJbL 

YLi  +  <Y  - 


-  1)M‘ 


1)M 


« 

72. 


dM 


M 


which  becomes  on  integration 


bL-~h-- h'W1* 


K 

1  +  (Y  -  l)M^/2j 

mJ 

1  +  <Y  -  UM3/2J 

YM2  YM3 

For  small  Mach  numbers,  this  can  be  written  as  (see  equations  (7)) 

,2 


bL 


1  *  (P^pp 


-  (127!'i  ln(P4/Pi)' 


(11) 


(1?) 


2 

This  equation  differs  from  (10a)  by  the  {3  coefficient,  which  ia  due  to  the 

area  changes,  and  by  the  logarithm  term.  Again,  this  term  is  small  com- 
2 

pared  to  the  p  term,  except  at  small  pressure  ratios.  Henceforth,  we  refer 
to  equation  (12),  without  the  logarithm  term,  as  Darcy's  equation. 

When  the  pu  du  term  is  retained,  we  obtain  in  place  of  equation  (11) 

(see  [4],  Chapter  6) 

bL  =  F(M2)-  f(m2)  ,  (13) 

where 


F(M2)  = 


yM 


h  + 

r  m2  i 

\  2y  T 

1  +  (y  -  i)M  /2. 

(14) 


which  is  shown  in  Fig.  7  for  y  =  (7/5).  For  comparative  purposes,  this  figure 

2  2 
also  shows  1/yM  ,  which  closely  approximates  F  for  small  M  .  For  small 

Mach  numbers,  equation  (13)  becomes 


*.T  a2  1  "  (p4/pl)  ,  /y  +  1\,  ,  ,  ,2 

bL  =  P  - - +  {%-)  Mp4/Pl) 

yMj  1 


(15) 


which  again  differs  from  Darcy's  equation  by  the  logarithm  term.  This 
difference  is  also  negligible  except  at  small  pressure  ratios.  It  is  interest¬ 
ing  to  note  that  equations  (12)  and  (15)  are  identical  to  the  isothermal  results, 
equations  (10b),  when  y  =  1. 

2 

Darcy' 8  equation  shows  that  p  can  be  absorbed  into  b  in  the  incom- 

2 

prescible  case.  In  this  case,  only  one  constant,  b/p  ,  characterizes  the 


gas- «  ©lid  interaction.  This  is  not  surprising,  since  [5]  develops  a  theory  in 
which  (5  de<sa  not  enter  explicitly.  When  compressibility  becomes  important 
this  simplification  is  no  longer  possible  and  both  b  and  (3  are  separately 
significant  parameters. 

We  have  shown  that  Darcy's  equation  can  be  considered  as  the  incom¬ 
pressible  limit  for  both  the  isothermal  and  adiabatic  cases,  providing  p^/p^ 
is  not  too  small.  An  estimate  of  a  lower  limit  for  p^/'p^,  denoted  by 
{p4/pt  )min>  can  be  obtained  by  assuming  the  magnitude  of  the  logarithm  term 
is,  say,  1/10  that  of  the  other  term.  This  results  in 


(p4/p,> 


min 


exp 


£ 


L  10(y+i)M“ 


(16) 


2  2 

where  the  (p^/p^)  quantity  in  the  6  term  has  been  neglected.  Fori  example, 

equation  (16)  results  in  (P4/Pi)min  =  0.044  when  y  =  (7/5),  p  =  (1/2),  and 

Mj  =  0.  001,  Hence,  if  0.  044  <  (P^/Pj).  Darcy's  equation  can  be  usbd  in  place 

of  equation  (15),  providing  all  Mach  numbers  are  small.  Equation  (16)  shows 

that  small  values  ofip/Mj  result  in  values  of  (P4/Pj)min  near  unity,  in  which 

case  Darcy's  equation  is  no  longer  useful. 

A  criterion  for  the  use  of  incompressible  flow  theory  can  now  be  derived. 

2 

In  doing  so,  we  approximate  F  by  1/yM  and  freely  use  equations  (7).  Since 
.he  largest  of  the  Mach  numbers  is  M^,  the  flow  is  considered  incomprecsible 

7 

if  F(M^)  2  10  or  M^  -  0,  239  (see  Fig.  7),  From  equation  (13),  we  have 


— j  =  -  bL  £  10 

yM.j  V^2 
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which,  after  some  manipulation.,  results  in 


• 


(P4/PA)2  *  {P^Ppfn 


The  flow  is  thus  approximately  incompressible  if  (P4/p^)  is  greater  than 
^/P^inc*  ^is  rcsu*t  should  not  be  confused  with  (16),  which  concerns 
the  relative  importance  of  the  logarithm  term  in  equation  (15). 

The  Compressible  Momentum  Equation 
In  the  remainder  of  this  paper,  we  deal  only  with  the  compressible, 
adiabatic  result,  equation  (13).  This  equation,  in  conjunction  with  equa¬ 


tions  (2),  shows  that 


bL  =  fQ(Mj,p,  y,p4/p1) 


An  important  conclusion  based  on  equations  (1)  and  (18)  is  that  for  given  value 

of  P,  y,  p4/p1,  and  p^Pj,  the  mass  flow  rate  p^u^  adjusts  to  conserve 

momentum.  Any  change  in  the  flow  rate  through  a  plate  thus  requires  a  new 

p4/p^.  For  fixed  values  of  p,  y,  and  p4/p^,  an  alternative  interpretation  is 

2 

that  momentum  conservation  relates  and  bL.  In  fact,  this  relation  is 
single  -valued  as  shown  in  the  Appendix. 

As  noted  earlier,  the  two  Fanno  curves  coincide  when  p  =  1.  We 


therefore  obtain  the  relation 


».  V.  p4/p,) «  f(mJ)  -  f(m|)  . 
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in  order  to  maintain  constant  conditions  at  1  and  4. 


2  _ 2 

We  conclude  this  section  by  examining  the  limit  -*  1,  with 
2  _  ^  2 

and  (PjM^/p4)  =  (p^M^/p^)  .  In  this  situation, 

bL  =f1(Mj/p2,Y)  =  F(M^)- F(l)  ,  (20). 

as  is  evident  from  equations  (2b),  (13),  and  (14).  This  result  is  shown  as 

Fig.  8  for  y  =  (7/5).  Hence,  for  given  values  of  y  (e.  g. ,  7/5),  p,  and  p^/p^  we 

_ 2 

can  determine  bL  such  that  =  1.  To  do  this,  first  determine  by 

Fig.  4,  then  (p^M^/p^)  by  Fig.  5,  and  finally  bL  by  Fig.  8. 

We  may,  of  course,  specify  y,  P>  and  bL  and  then  determine  a  value 

for  p^/pj.  denoted  in  accordance  with  our  notation  by  (p^/pj),  such  that 

=  1.  Figure  9  shows  (p^/p^)  versus  bL  for  y  =  (7/5)  and  various  values 
2 

of  p  ,  including  p  =  1.  Values  of  P^/Pj  greater  than  (P^TPj)  result  in  <  1; 
values  less  than  (p^/p^)  are  considered  in  the  next  section.  All  the  curves 
in  Fig.  9  are  parallel  because  (M^/p  )  =  (p^7p|)  g,  where  g  is  a  function 
only  of  p  and  y.  The  solid  curve  is  for  both  P  =  0  and  p  =  1.  As  P  increases, 
m  at  first  increases,  reaches  a  maximum  at  about  P  =  0.  6,  and  then 
decreases  to  the  P  =  1  value.  This  behavior  can  be  shown  to  be  due  to  the 
curvature  of  the  path  that  represents  the  3  to  4  process,  as  shown  in  Fig.  2. 
The  process  twice  intersects  some  constant  pressure  lines,  thereby  resulting 
a  double -valued  behavior  for  (p4/p^).  Since  the  curves  are  so  close  together, 
the  pressure  ratio  for  choking,  (p^/p^),  is  insensitive  to  p. 


■3* 


Finally,  for  large  bL,  i.  e. ,  for  a  thick  plats,  we  have  approximately 

i  fl2 

bL  ss - k  *  ■*  s 

ym|  ymJ 

2 

and  hence  ln(p^/pp  is  proportional  to  In  bL,  as  is  evident  from  Fig.  9. 

2  2 

Note  that  bL  * (P  / yM^)  is  essentially  Darcy's  equation  when  P4/pj  is  small. 
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CHOKE.D  FLOW:  M3  =  1. 

When  (p4/pj)  <  (p4/p1),  equation*  (,2c),  (3),  (4)>  .aijS  (5) .no  longer  apply, 
and  state  4  cannot  be  directly  related  to  state  3  a*  before.  These  equations 
are  replaced  by  =  i,  and  condition*  at  i,  2,  and  3  ape  unaffected  by 
changes  in  p^.  Since  conditions  at  1  and  2,  as  determined  by  equations  (2b) 
and  (20),  are  fixed,  we  concentrate  on  conditions  at  4,  where  is 
determined  by  equation  (2a). 

The  analysis  is  again  conveniently  illustrated  by  a  Mollier  diagram 
(Fig.  10).  State  points  1,  2,  and  3  are  the  same  as  in  Fig.  2,  except  that 
3  is  at  the  "nose"  of  the  Fanno  curve.  The  limit  “*.i,  is  shown  as  4.  For 
pressure  ratios  slightly  less  than  (p^/pj).  state  4  is  to  the  right  of  4  on  the 
Fanno  curve  through  state  1.  State  4  thus  moves  in  a  continuous  manner 
along  the  Fanno  curve  as  P^/pj  decreases. 

1  If  p^  is  further  decreased,  state  4  is  at  the  "nose"  of  the  curve  where 

K  jjg 

MJ  =  1,  and  the  corresponding  pressure  ratio  is  designated  by  (P^Pj)  •  A 
furcher  decrease  in  p^/p^  results  in  supersonic  flow,  since  >  1.  For  a 
sufficiently  small  pressure  ratio,  denoted  by  (p^/pj),  conditions  at  4  are 
realhed  by  an  isentropic  expansion  from  3  to  4,  as  shown  in  Fig.  10. 

Solutions  do  not  exist  for  pressure  ratios  smaller  than  (p^/pj),  which  is 

1  2 

referred  to  as  the  limiting  pressure  ratio.  This  limit  results  from  our 

2  Tl  !  TT" 

This  limiting  pressure  ratio  probably  cannot  be  attained  in  practice  due  to 
the  effects  of  turbulence  downstream  of  the  plate.  A  smaller  downstream 
reservoir  pressure  than  the  limiting  one  should  result  in  a  flow  at  the  end  of 
the  duet  analogous  to  the  flow  external  to  an  underexpanded  nozzle. 
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insistance  on  &  constant-area  duct.  A  snoaUer  pressure  ratio*  for  example, 

.  ^requires  an  expanding  flow'  and  an  increasing  duct  cross-sectional  area 
downstream  of  the  plate* 

it  should  be  realised  that  the  flow  may  be  quite. turbulent  downstream 
of  the  plated  In  this  situation,  state  4  does  net  oe-cur  until  a  uniform  duct 
flow  is  achieved.  Thus,  state  4  is  frequently  downstream  of  the  plate  rather 
than  at  its  surface. 

Conditions  at  4,  when 

•I 

I  (P4/Pt)  S  (P4/Pt)  SCpi/Pj) 

are  determined  by  combining  equations  (2a>,  (2b),  and  (13)  to  yield 

F{hi|)  sfcL  +  F(i)  ,  (21a) 

m4 =rr^7)  |[*  +  2(v  ■  1)(Pi/p4^2  gi(m2)]1/z  -  *J  •  {21b) 

Thus,  M4  depends  only  on  y,  bL,  and  (Pj0/pt4)-  The  pressure  ratio  (p^Pj)*, 
as  a  function  of  bL,  y,  and  P,  is  obtained  by  setting  M4  =  1.  Figure  9  shows 
[(P4/Pj)  1  versus  bL  for  y  =  (7/5)  and  various  0  ,  including  0  =  1.  As  with 
(i^7),  these  curves  are  also  parallel.  One  difference,  however,  is  that 
decreasing  0  results  in  a  smaller  (p4/pj)  for  a  fixed  bL.  This  is  evident 
once  it  is  realized  that  decreasing  0  separates  the  Fanno  curves,  while  a 
constant  bL  means  that  the  Fanno  curve  through  2  is  kept  fixed. 
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The  Mach  number  at  the  limiting  pressure  ratio  is  determined  by 
the  isentropic  relation 

=  [2/(y  +  i>]<Y+1)/(V‘1>  P2  . 

where  equations  (21)  still  apply.  These  may  be  solved  for  (p^/Pj,)*  which  is 

2  ’ 

also  shown  in  Fig.  9  for  y  =  (7/5)  and  various  {3  ,  including  P  =  1,  The 
variation  of  (p^/pj)  with  P  is  much  greater  than  for  the  other  pressure  ratios 
shown  in  Fig.  9.  We  also  see  that  all  the  curves  are  again  parallel. 
Basically,  this  occurs  because  the  same  momentum  equation  (13)  and  Fanno 
curve  equation  (2a)  apply  to  all  the  curves. 


DISCUSSION 


Three  principal  topic*  are  discussed  in  this  section,  all  of  which  bear 
on  experimental  verification.  First,  the  various  flow  regimes  are  described 
in  terms  of  the  pressure  ratio.  An  appropriate  generalization  of  the  volumetric 
flow  rate  versus  pressure  differential  plot  is  then  given.  The  third  topic  is 
a  discussion  of  th©  physical  implications  of  the  drag  coefficient  b. 

/ 

There  are  four  flow  regimes  that  can  be  characterised  by  the  pressure 
ratio  p^/pj.  The  first  is 

(p4/pl>ii.c  s  <p4/pl>  s  1  ’ 

where  (P^Pj)^  1®  given  by  equation  (17).  This  is  the  incompressible  regime 
where  Darcy's  equation  is  usually  valid.  The  next  regime 


(P4/Pt)  S  (P^Pj)  <  <P4/Pi) 


xnc 


is  characterised  by  compressible  subsonic  flow  everywhere.  The  third 
regime 

(P4/P1)*  5  (p4  ^Pj )  K  (P^Pp 


requires  a  choked  flow,  with  subsonic  flow  at  4.  The  final  regime 

$ 

(P4/P^  s  (p4  /pA )  <  (P4/Pj) 

also  requires  a  choked  flow,  with  supersonic  flow  downstream  of  the  plate. 

Of  the  four  regimes,  there  can  be  no  doubt  about  the  characteristics  of  the  first. 
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i 

* 

The  other*  *tiU  require  experimental  verification.  Thi*  is-  particularly 
true  of  the  last,  since  supersonic  flow  has  not,  as  yet,  been  associated 
with  flow  through  a  porous  medium.  This  regime  may  be  the  easiest  to 
verify,  however,  since  standard  techniques,  such  as  the  Schliere'n. optical 
method,  can  be  used  to  detect  the  presence  of  supersonic  flow  downstream 
of  the  plate. 

Incompressible  flow  through  a  porous  medium  is  frequently  represented 
by  a  linear  relation  between  the  volumetric  flow  rate  and  the  pressure  dif¬ 
ferential.  This  representation  does  not  hold  for  compressible  flow.  More 
appropriate  would  be  a  plot  of  M2  versus  (p^/pj)2,  as  shown  in  Fig.  11. 

(Once  Mj  and  (P4/Pj)  are  known,  is  readily  determined  by  equation  (2a). ) 
This  figure  is  for  02  =  0.  2,  y  =  (7/5),  mid  various  bL.  The  solid  curves 
are  for  compressible  flow,  where  the  (A)  denotes  the  pressure  ratio  , 
(P4/Pj)jnc»  an^  the  (o)  denotes  the  onset  of  choking.  These  curves  terminate 
at  (p^Pj)*  The  dashed  line  is  the  incompressible  result,  equation  (15),  while 
the  doubly  dashed  line  is  Darcy's  equation,  which  differs  from  (15)  when 

(P4/Pj)  5  ^4^1  ^min'  This  difference  shifts  to  smaller  pressure  ratios  as 
bL  increases  in  accordance  with  (16). 

For  bL  equal  to  1  and  10,  we  see  that  the  compressible  and  incompres¬ 
sible  results  depart  at  (P4/P1)inc-  For  bL  equal  to  102  and  (p4/pt)2  <  0. 1, 

z 

only  small  changes  in  Mj  are  necessary  for  large  changes  in  and  P4/pj. 

Since  this  is  true  for  both  the  incompressible  and  impressible  results,  the 

two  remain  quite  close.  For  a  thick  plate,  i.  e. ,  large  bL,  Darcy's  equation 

2  2 

adequately  predicts  versus  (P4/Pj)  ,  as  noted  earlier.  The  incompressible 
solution  does  not,  of  course,  correctly  predict  M4  or  (p4/pj). 

Z  - 


•  Experimental  data  will  not  necessarily  follow  a  single  bL  curve  *s  flow 
conditions  vary.  Deviations  occur  because  b  generally  varies  with  flow 
conditions.  It  is  therefore  necessary  to  explore  the  relationship  of  b  to  the 
other' parameters  of  the  problem*  since  experimental  verification  depends  to 
some  extent  on  this  relationship. 

Any  experimental  investigation  must  deal  with  at  least  three  types  of 
parameters.  The  first  are  those  that  characterize  the  gas,  such  as  y  and 
the  viscosity.  The  second  are  those  that  characterize  the  porous  medium, 
such  as  p.  Finally,  there  are  parameters  associated  with  the  specific  flow 
situation,  such  as  PjU^.  All  parameters  used  so  far,  except  b,  have  a 
straightforward  physical  interpretation  and  are  experimentally  determinable. 
Part  of  the  difficulty  with  b  is  that  it  depends  on  the  gas,  porous  medium, 
and  flow  conditions.  This  is  evident  from  the  relation  b  =  (Zrj/puBg)  given 
earlier. 

The  drag  coefficient  was  introduced  with  the  expectation  that  it  is 
independent  of  p,  which  must  be  specified  as  well.  In  addition,  it  is  reasonable 
to  expect  that  b  does  not  depend  on  the  molecular  weight  of  the  gas  or  on  y. 

The  first  exclusion  stems  from  dimensional  analysis,  while  the  second  is  due 
to  the  explicit  appearance  of  y  in  the  theory.  In  order  to  define  more  clearly 
the  role  of  b,  we  proceed  by  comparing  the  flow  in  a  porous  medium  with  that 
in  a  pipe  with  roughened  walls.  This  type  of  comparison  is  subject  to 
criticism,  as  [i]  points  out,  and  our  conclusions  are,  at  best,  heuristic. 

We  nevertheless  briefly  pursue  this  course  in  order  to  introduce  concepts 
that  may  be  experimentally  useful. 
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According  to  this  analogy,  b  should  depend  on  some  average  pore 
diameter  d  (or  equivalently  on  the  specific  surface),  and  on  the  Reynolds 
number  =  {p^u ^d/rj^ ).  Thus,  b  =  (Zd/B^R^)  and  if  d/B^  is  constant,  b  is 
inversely  proportional  to  the  Reynolds  number,  a  result  typical  of  laminar 
flows.  For  sufficiently  large  values  of  R^,  however,  due  to  turbulence, 
d/BQ  is  not  necessarily  constant.  The  drag  coefficient  may  also  depend  on 
the  Mach  number,  although  this  dependence  is  believed  to  be  weak  {see 
[4],  Chapter  6).  For  a  specific  medium,  b  may  also  change  with  time  due 
to  surface  changes  such  as  adsorption.  In  addition,  b  may  depend  strongly 

on  certain  characteristics  of  t^e  material,  such  as  qhe  relative  roughness 

\ 

(see  [4],  Chapter  6).  For  sufficiently  large  flow  rates,  the  flow  is  not  only 
turbulent,  but  b  may  be  nearly  independent  of  R^.  In  this  circumstance,  b 
depends  primarily  on  the  relative  roughness  alone.  This  regime  is  ideal  for 
experimentally  verifying  this  analysis,  providing  it  occurs  before  the  flow 
chokes. 

Even  if  b  does  not  have  precisely  the  behavior  described  above,  experi¬ 
mental  verification  is  still  possible.  For  example,  use  can  be  made  of  the 

fact  that  once  choking  has  occurred,  the  flow  from  2  to  3  does  not  change 

3 

and  bL  remains  constant.  Furthermore,  in  the  compressible  regime  when 
<  1,  problems  associated  with  b  might  be  avoided  by  using  different  plate 
thicknesses  of  the  same  porous  material.  For  engineering  purposes,  plots  of  b 
versus  the  pertinent  parameters,  such  as  R^,  are  of  course  necessary. 

3 

We  are,  of  course,  allowing  only  p^  to  vary. 
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In  conclusion,  note  that  although  we  have  dealt  exclusively  with  the 
pressure  and  Mach  number,  other  quantities  are  easily  determined.  For 
example,  the  temperature  at  4  is  known  for  all  flow  situations  from 


9 


since  the  stagnation  temperature,  approximately  given  by  Tj,  is  a  constant. 
Another  useful  quantity  is  the  stagnation  pressure  p-.  At  1  this  is  approxi¬ 
mately  equal  to  p^,  but  at  4  it  is 


£04 

?4 


/'  « 

=  [i  +>r^i)M*/2] 


iV/(v-l) 


I 


and  consequently,  (p^/p^)  <  (P^/Pot^  Tlli®  difference  in  static  and  stagna¬ 
tion  pressure  ratios  is  not  necessarily  small,  e.  g. ,  with  y  =  (7/5)  and 
M4  =  1  we  have  (P04/P01)  =  1. 892  (p4/pt ). 
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APPENDIX 


If  (sbL/SMj)  it  non-zero*  then  bL  and  are  uniquely  related.  This 


derivative  it  [  tee  equation  (18)] 


where  the  derivative!,  such  at  (sM^/SMj),  are  determined  by  equations  (2) 

% 

and  (4).  We  thus  obtain 


8(bL)  1  [  1  _  JL\ 

8M^  yMj  \m|  Mg/ 


where 

.  (i-u\)(i+yml)(i-ul) 

3  [i  +  (V  -  1  -  VPJMgJ  (l  +  V  M^)  Tl  +  (y  -  i)M*J 

2  2 

and  Mg  and  M^  are  related  by  equation  (4).  It  is  easy  to  see  that  H(0)  =  i  and 
H(l)  =  0.  One  can  also  show  that  |dH/dMgj  £  0  for  0  £  Mg  £  '  Consequently, 
H  £  1  and  since  M2  <  Mg,  we  have  ^8bL/8M^j  <  0,  and  bL  and  M^  are  uniquely 
related.  Furthermore,  since  |8bL/8M^  j  is  negative,  increasing  M^  means 
bL  must  decrease  if  p^/pj  is  to  remain  constant.  In  other  words,  we  get 
the  result  that  if  the  pressure  ratio  p^/p^  is  to  remain  constant  as  the  mass 
flow  rate  is  increased,  the  dimensionless  product  (bL)  must  decrease. 
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POROUS  PLATE 
CROSS  SECTION); 


(P4/Pl>2 


Fig.  6.  Comparison  of  the  Terms  in  Equation  (10b); 
Y  =  7/5,  M?  =  0.  01 
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